Long-range effects in granular avalanching 
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O^l ■ We introduce a model for granular flow in a one-dimensional rice pile that incorporates rolling 

j—t ' effects through a long-range rolling probability for the individual rice grains proportional to r~ p , 

r being the distance traveled by a grain in a single topling event. The exponent p controls the 
\^ • average rolling distance. We have shown that the crossover from power law to stretched exponential 

behaviors observed experimentally in the granular dynamics of rice piles can be well described as a 
long-range effect resulting from a change in the transport properties of individual grains. We showed 
that stretched exponential avalanche distributions can be associated with a long-range regime for 
f < p < 2 where the average rolling distance grows as a power law with the system size, while 
power law distributions are associated with a short range regime for p > 2, where the average rolling 
C/} ' distance is independent of the system size. 
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PACS numbers: 64.60.Lx, 46.10.+Z, 05.40. +j 



There is nowadays a massive evidence of scale invariant phenomena in nature. They appear in systems as diverse 
G ■ as geological (e.g., earthquakes), climatic (e.g., atmospheric turbulence), granular flow (e.g., rice piles, the topic of 
i-^j | this work) and biological (e.g., biological evolution, cell growth, among many others). In many of these cases no 
particular tunning is perceived. In 1987 Per Bak and colaborators JIJ] advanced the hypothesis that this is so because 
the microscopic dynamics of the system makes it to spontaneously evolve towards a critical, scale-invariant state. This 
i^i - is known today as self-organized criticality (SOC). To illustrate the idea they introduced a sand pile model pj| that 
quickly became the paradigm ofSOC models. In this model a pile is gradually built by adding individual grains one 
■ by one into an open system. After a transient period the pile reaches a stationary critical state, where the slope of the 
pile fluctuates around a stationary value. At this stage avalanches of all possible sizes are observed, giving a power 
^ law avalanche size distribution. This is the SOC state. The experimental search for power laws in the avalanche size 
distribution in real granular systems presented a great challenge Some years ago the Oslo grou p @ ,[l0| found 

evidence of SOC in controlled experiments on the granular dynamics of rice piles. In 1996 Frette et al [0 performed 
an experiment where elongated rice grains were added slowly in a narrow gap between two plates. They found that 
the avalanche size distribution for grains with large aspect ratio presents a power law behavior, while a stretched 
exponential behavior is observed for rounder grains. These experimental results showed that SOC is not insensitive 
to the details of the system such as the shape of the grains and that inertia effects play an important role in the 
relaxation dynamics. The Oslo group introduced a "rice-pile" model to phenomenologically describe their experiments 
I 1 [jnjjl' This model describes the phenomenology of rice piles presenting SOC behavior. Moreover, further studies 
showed that it can be exactly mapped to a model for interface depinning jljj] . However, the Oslo model does not 
£h , present a stretched exponential distribution in any region of its parameter space. 

In a recent communication Head and Rodgers ]T^] presented a clustered based model for granular flow, which 
exhibits both stretched exponential and power law behavior over its parameter range. Through this model they 
showed that power law behavior appears to be related to the coherent motion of clusters of grains that slide along 
the surface, while the stretched exponential behavior results from a fluid like motion associated with the rolling of 



isotropic grains. However, the motion of individual grains is implicitly incorporated in this model, where the basic 
dynamical variables are related to cluster properties. Hence, it is interesting to investigate in a more microscopic 
level how the crossover from power law to stretched exponential behavior emerges as a consequence of a change in 
the transport properties of individual grains. 

In this work we present a model where the ability of individual grains to roll a distance r is described by a long-range 
rolling probability of the form 

p(r)= U/r<M/ l<r<A 
10 otherwise 

where A is a normalization constant, N is the system size, r = 1,2, ...,N and < p < oo. The parameter p 
controls the average rolling distance which is expected to depend on the aspect ratio of the grains. The p — ► oo 
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limit corresponds to a nearest neighbors movement (P(r) = 6 r ,i), thus describing the case where the grains do not 
roll. In the opposite limit p = the grains can move to any site of the pile (or even drop off in a single jump) with 
equal probability. Hence, this limit can be associated to the idealized situation of perfect spherical grains without 
friction. Moreover, studies on different dynamical systems with similar long-ranged interactions [jl3]-[l5| have shown 
that, except from some possible rescaling factor, the behavior of almost all the relevant properties of the systems in 
the full range < p < 1 reproduce the corresponding ones for p = 0. This appears to be related to the behavior of the 
first moment of the distribution (Q), which in the present context represents the mean rolling distance of the grains. 
For large values of N we have that 
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We see that, for < p < 1, (r) is of the order of the system size, which means that a very large number of grains will 
drop off of the system in a single jump, for any system size. Since this behavior appears to be highly unrealistic we 
will concentrate our study on the p > 1 regime . 

Our model is based on the Oslo one [j§|-[l0[ and it is defined as follows. We consider a one dimensional lattice of size 
N (1 < i < N), each site i having associated an integer variable h(i) representing the local height of the pile. The 
local slope is then given by a(i) = h(i) — h(i + 1). The grains enter into the system from the left (i = 1) and may 
drop off at the rightmost site i = N + 1, imposing h(N + 1) = for all times. Every time the local slope cr(i) of a site 
i exceeds a local critical value a(i) > a c (i) the topmost grain at site i rolls r sites to the right with probability P(r) 
given by Eq.(|l]). Then, the heights of the sites i and i + r are recalculated as h(i) — > h(i) — 1 and h(i + r) — > h(i + r) + 1 
and the corresponding local slopes are modified accordingly. Each time a grain leaves a column i we assign ita new 
critical slope a c (i) = 1 or a c (i) = 2 with equalprobability. This process is repeated until all the local slopes satisfy 
o{i) < a c {i)- An avalanche starts when cr(l) > er c (l) and when it stops (i.e., < cr c {i) Vz) new grains are added 
until a new avalanche is initiated. Notice that in the limit p — > oo the present model reduces to Oslo model. 

Once the system attains the stationary state, fluctuations in the slope of the pile appear in the form of avalanches 
that redistribute the grains, thereby changing the profile. By considering the profile just before an avalanche is 
initiated and the final profile when it stops, the avalanche size is defined as the total energy dissipated between both 
profiles §. The local energy change in a site i is calculated as the difference between the potential energies in both 
profiles, where the units were chosen such that mg — 1. We study for different values of p the probability density 
P(E,N), P(E, N)dE being the probability that an avalanche with energy dissipation between E and E + dE will 
occur in a system of size N. Numerical data were smoothed using a local average procedure in order to diminish large 
statistical fluctuations present for high values of the ratio E/N. We verified for a wide range of values of p (even 
in the non physical region p < 1) that P(E,N) displays the finite-size scaling behavior P(E,N) = N^ 1 f(E/N), in 
agreement with the experimental data 0, as shown by the data collapse of Fig. [j] for p = 1.4 and p — 4. The scaling 
function f(x) depends on the value of p. We found two clearly distinct regimes. For p > 2 (hereafter referred as the 
"sliding regime") f(x) is almost constant for small values of x and presents a power law dependency f{x) ~ x~ a for 
large values of x, as shown in the example of Fig.|| for p = 4. The exponent a increases smoothly from a « 1 in the 
limit p — > 2 + to a = 1.33 ± 0.05 for values of p > 4, where it becomes independent of p. This value is consistent with 
that obtained by Head and Rodgers |0| for the cluster based rice pile model. 
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FIG. 2. Log-log plot of P(E, N)N vs. the energy dissipation E/N for TV = 1000 and p = 4. The same qualitative shape is 
observed in the whole regime p > 2. The straight line indicates the best fitting in the power law region fix) ~ x _a . 



For 1 < p < 2 (hereafter referred as the "rolling regime") the best fitting of the numerical data is obtained by a 
function of the form: 



f( x ) K ~z exp 



(3) 



as shown in the example of FigJ^ for p = 1.2. The parameters 7, a and xo depend on p, as shown in table Q, 
although the small variation of 7 in this regime suggests that it may not depend on p. Notice that this function seems 
to describe the behavior of the whole distribution, not just for large values of E/N (where it is dominated by the 
stretched exponential) but also for small values of it. 
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FIG. 3. Log-log plot of P(E,N)N vs. the energy dissipation E/N for N = 1000 and p = 1.2. The same qualitative behavior 
is observed in the whole regime 1 < p < 2. The dahed line corresponds to the best fitting using a function of the form 
f(x) oc x~ a exp(— (x/xo) J ) (see table | for the values of the fitting parameters). 



We have shown that the crossover from power law to stretched exponential behavior observed in the granular 
dynamics of rice piles can be well described as a long-range effect resulting from a change in the transport properties 
of individual grains. Within this scenario two distinct regimes appear, regarding the qualitative behavior of the 
avalanche distribution: a short-range "sliding" regime and a long-range "rolling" one. A direct inspection of snapshots 
of the profiles obtained during the simulations shows that, indeed grains group into clusters that move coherently 
in the sliding regime, while in the rolling regime they move almost independently. These behaviors are consistent 
with what is observed experimentally and with Head and Rodgers results Details of these simulations will 
be published elsewhere. This change in the transport properties of individual grains is related to different scaling 
properties of the average rolling distance with the system size: while in the rolling regime d scales as a power law 
d ~ N 2 ~ p , in the sliding regime d does not depend on N . 

In the whole sliding regime p > 2 the system displays the qualitative behavior expected for a short-range model, 
that is, a power law for large values of the energy density E/N and constant for small values of it. However, it is 
worth noting that for intermediate values of p (2 < p < 4) the characteristic exponent a of the power law distributions 
depends on p. while for large values of p it becomes independent of p. This results suggests that non-universal behavior 
may be expected for intermediate values of the grains aspect ratio. The value p — 4 is only indicative and it seems 
not to be a critical value of p separating the universal from the non-universal behaviors. 

In the rolling regime the whole avalanche distribution is well described by a stretched exponential modulated by a 
power law. This correction to the usual stretched exponential behavior is only important in the small energy density 
region of the spectrum, since for x 3> xo the function is always dominated by the stretched exponential. Notice that 
for 1.4 < p < 2 the exponent a changes of sign taking small absolute values, suggesting logarithmic rather than power 
law corrections to the stretched exponential in this region. Indeed, a slightly better fitting is obtained with a function 
of the form 

f(x)ocln(x)exp(-(j-y^ (4) 

without any significant change in the fitting parameter values of xq and 7 shown in table [|. 

From table | we see that the best agreement |?J with the experimental values xq — 0.45 ± 0.09 and 7 = 0.43 ± 0.03 is 
obtained for p = 1.01. Both functions Eq.(^) and cease to give good fittings of the simulation data as we approach 
the border values p — ► 1 + and p — ► 2~ respectively. This fact, together with an increase in the statistical fluctuations 
in the neighborhood of those values, makes very difficult to determine how the crossover from stretched exponential 
and power law occurs around p = 2. Possibly the calculation of the distributions of other quantities (like transit 
times) may be useful. Some work along this line is in progress. 
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0.15 ±0.05 


0.23 ±0.05 


0.35 ±0.05 


1.2 


0.15 ±0.05 


0.20 ±0.05 


0.35 ±0.05 


1.4 


0.05 ±0.05 


0.15 ±0.05 


0.33 ±0.05 


1.6 


-0.03 ±0.05 


0.11 ±0.05 


0.31 ±0.05 


1.8 


-0.09 ±0.05 


0.06 ± 0.05 


0.28 ±0.05 



TABLE I. Fitting parameters using the function (j3|) in the rolling regime 1 < p < 2 
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